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Abstract 

It is known that many peculiar nonlinear vibration problems in impacting systems are caused by 
grazing incidences. Such bifurcation phenomena are normally investigated through the Poincare map. 
The discrete-time map of a simple impact oscillator was derived by Nordmark, which showed that there 
should be a square-root singularity in the Jacobian matrix close to the grazing condition. In this paper 
we show that the square root singularity will be expressed only in the trace of the Jacobian matrix, while 
the determinant remains invariant across the grazing condition. 
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1 Introduction 

Mechanical systems with impacts between elements occur frequently in engineering practice. Vibration 
problems in such systems essentially hinge on the dynamics of a moving body, possibly in a mass-spring- 
damper combination, impacting with a rigid stop. It is known that much of the dynamical phenomena in 
such systems stem from the conditions pertaining to the situation when one body just grazes the other. That 
is why much attention has been given to the grazing condition, and the bifurcation phenomena resulting 
from that [HEiEiii]. 

It is convenient to analyze bifurcation phenomena in any dynamical system by obtaining a discrete time 
model or map by the method of Poincare surface of section. The structure of the obtained map determines 
the dynamics of the physical system. Thus, in understanding the dynamics of impacting systems, researchers 
have tried to obtain the structure of the map. It is obvious that the map should be piecewise smooth, since 
the map for non-impacting condition and that for the impacting condition should be different, and the two 
should be separated by the grazing condition. Most important in the respect is the question: What is the 
structure of the map in the neighborhood of a grazing orbit? 

Nordmark [4] first addressed this question. He showed that in the non-impacting side the map is linear, 
while in the impacting side it has a square-root term. This implies that the derivative of the map approaches 
infinity as the grazing condition is approached from the impacting side. This results in an infinite stretching 
of the state space—which has come to be known as “square-root singularity.” Many researchers studied the 
behavior of the impact oscillator based on this map [SJ1SJ171IH]. 

In a two-dimensional oscillator (mass-spring-damper type), a sampling in synchronism with the external 
forcing finction yields a two-dimensional map. Nordmark derived the condition on the whole Jacobian matrix, 
but not on the elements of the matrix. In the meantime, many other switching dynamical systems—most 
notably the power electronic circuits—were found to exhibit a new class of bifurcation that occurs when a 
fixed point crosses the border between two smooth regions in a piecewise smooth map. The development of 
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the theory of such border collision bifurcation based on a normal form . The normal form is expressed as 



( 1 ) 


where is the trace and <5^ is the determinant of the Jacobian matrix Jl of the system at a fixed point in 
one side and is the trace and 5r is the determinant of the Jacobian matrix J_r of the system evaluated 
at a fixed point in the other side. 

The natural question in relation with the impact oscillator was: How does the trace and the determinant 
change as an impact oscillator is driven from a non-impacting state to an impacting state with the change 
of a parameter? 

In this paper we probe this issue, and analytically prove that the determinant must be invariant while 
the trace alone should exhibit the square-root singularity. 


2 Impacting Hybrid System Description 


An impacting system (Fig. [T]) is governed by a set of ordinary differential equations (ODEs) coupled with a 
set of reset maps as 


i = F(x), if X € 5+ 

xH>i?(x), if xGE (2) 

where, = {x : i/(x) > 0} and E = {x : i/(x) = 0}. i?(x) is a smooth function, zero set of which defines 
the hard boundary E. The flow given by ([2]) is restricted only in the region S'^ U E. 



Figure 1: Upon impact, the velocity instantly reverses while the position remains the same. Thus the state 
instantly jumps to a new position. 


Let us now define the normal velocity u(x) as the rate at which the trajectory approaches the impact 
boundary. It is given by 


dH dH dx 









Similarly the normal acceleration a(x) of the flow with respect to the boundary is 


a(x) := {H^F)^F. 

We may now be more specific about the form the reset map i?(x) takes. To that end, we observe that the 
reset map has to be a smooth function of the normal velocity v(x) and furthermore i? maps to itself when 
grazing occurs. Since at grazing the normal velocity with respect to the boundary becomes zero (u(x) = 0), 
the reset map can be formulated as 

i?(x) = X + W(x)u(x) (3) 

where W is a smooth 2x1 matrix. 


3 Grazing and Discontinuity Mapping 

Grazing occurs when a trajectory becomes tangent to the discontinuity boundary S, as shown in Fig. [21 A 
point X = X* is called a regular grazing point if it satisfies the conditions 


H{x*) = 0 

z;(x*) = 0 
a(x*) = a* > 0 

In addition the scalar function H{x) is assumed to be well defined at x = x*, i.e., i7x(x*) ^ 0. 



Figure 2: Grazing of discontinuity boundary. 


For the part of the flow that does not have any impact with the discontinuity boundary the mapping is 
given by the ODE only. Whenever there is an impact with the boundary, the reset map comes into action and 
there is a discontinuity in the flow. The discontinuity near grazing is of particular interest. Special kinds 
of mapping have been proposed to account for this discontinuity [10]. In this present case the zero-time 
discontinuity mapping (ZDM) is dealt with. 

Let us consider the situation as shown in Fig.|3l There is an orbit which grazes the discontinuity boundary E 
at a point x* at some point of time to- Now let there be another trajectory (X 0 X 1 X 2 X 4 ) close to the grazing 
orbit. Let us back-trace the trajectory, governed by the ODE as in from the point X 2 to the point X 3 
such that the time taken by the trajectory to reach from xq to X 2 is the same as would have been taken by 





the flow to reach from X 3 to X 2 . Thus, we can consider the systems’ dynamics as if the switching boundary 
were not there. In that case we have to assume an instantaneous jump of the state from xq to X 3 . The ZDM 
is defined as the mapping xq X 3 . 


u 


k 

I 



(o;o) 


Figure 3: Zero-time discontinuity mapping near grazing. 


Let (/)(xo, t) be the flow obtained as the solution of the ODE starting from xq, i.e. (/>(xo, 0) = xq. It has 
been shown in m that the form of the ZDM, excluding higher order terms, is 

xs = xo - W* {V2^)y (4) 

where, W* = IT(x*), and y = y'—iLmin(xo). iLmin(xo) is defined as the minimum value of with 

the smallest |t|, i.e., the lowest point that the trajectory would have reached if the switching boundary were 
not there. Obviously for the situation as described in Fig. |3l iLmin(xo) will be negative except for the case 
when Xq is the same as x*. 

Now, let us consider a periodic orbit which has an intersection with the discontinuity boundary very close 
to the grazing orbit, as shown in Fig. |4l The stroboscopic Poincare map in this case is Pg = P 20 ZDM 0 P 1 , 
where Pi is the map that takes a point on the Poincare plane and maps it to the discontinuity boundary S 
by evolution through the ODE in ([2]), and P 2 is the map that takes a point on the discontinuity boundary 
E and maps it back to the Poincare plane via the ODE. The form of this stroboscopic map can be derived, 
in first order approximation, as: 


Pi 

ZDMoPi 


X !->• A^ix 

X NiX - v^v'-FImin(iVix)VF* 

X ^ N2N1X - V^V-Hmin{Nix)N2W* 


P2OZDM0P1 


( 5 ) 







where, Ni := ^|x=xo, ^2 := ^|x=xo- 

On the Poincare plane, we define Xq as the origin, and the stroboscopic map in this case is defined as the 
map which takes the the initial deviation on the Poincare plane, (x — Xq), and maps this deviation again on 
the Poincare plane. Since Xq maps to the point N 2 N 1 X.Q, this map takes the form 

{x-Xo) 1 -^ N2Ni{yL-Xo) - V^y^^Il~JN^N2W* ( 6 ) 

Also iJinin(A^ix) is linearized about grazing to the form 

i?min(Aflx) = i?xiVl(x - Xo) + H.O.T 



Figure 4: A grazing periodic orbit and stroboscopic map. 


As can be seen from the expression of the stroboscopic map, and also of ZDM alone, a square-root term 
V—.ffmin is present which accounts for the square-root singularity when the Jacobian of the map is considered. 
The next section deals with this Jacobian and the square-root singularity therein. 


4 Investigating the Trace and the Determinant of the Jacobian 
for the Singularity 


The Jacobian of the stroboscopic map near grazing would be 

J = N 20 JZDM 0 N 1 (7) 


where Jzdm is the Jacobian of the ZDM given by 


7 T , ATT W*H-x_ 

JZDM = ^=I+V^.^^— 


( 8 ) 


To arrive at the particular forms W* and iJx would take, let us concentrate on the one degree-of-freedom 
impact oscillator (shown in figure 5). The mass (assumed to be unity without any loass of generality) is tied 
with a spring-damper element and is acted upon by an external force g{t). At a distance a from the mass 
there is an impacting wall so that for u < a the motion of the mass is governed by second-order differential 
equation 

(Pu 


du 










and at M = (T the reset map R is applied. This system is a two-dimensional system, i.e. x S Considering 
the velocity of motion of the mass as v = the state vector can be written as x = {u,v)'^. 



a 

Figure 5: A one degree-of-freedom impact oscillator. 


The equation of the discontinuity boundary S in the present case is 


H (x) = H{u,v) = a — u. 


Thus we have ^ = 0, and hence, 

ifx = (/ii 0). 


( 10 ) 


where hi = 

The reset map is i? : E i—>■ E, where i?(x) has the form given in ([3]). Note that the position u does not 
vary during the impact, i.e., the position of the mass just before the impact, it“, is same as that just after 
the impact, u'^, while the velocity of motion v changes. Thus from ([3|) 


VF = (0 W2f-, 


( 11 ) 


where W 2 is the constant of restitution. 
Using m and (HH) in (|8]) we obtain 


JzDM 


f 

1 0 
0 1 

1 0 
0 1 

1 0 


W*H^ 



( ") 

2V~tlniin 

\W2 J 

1 / 

' 0 

2\/~ttniin \ 

^ W2hi 


{hi 0) 


where 


W2hiy/2a* 


( 12 ) 


4.1 Investigating the Determinant for Singnlarity 

From dzD, the determiant of the normal form map near grazing is 


1^1 


|1V2||Tzdm||a^i| 

|iv 2 ||a^i| 





















since \ Jzdm\ = 1, from (IT^ . 

Since the singularity is only in the ZDM, and not in the maps Ni and A^ 2 , we conclude that the deter¬ 
minant of the normal form map does not contain the square-root singularity, and remains invariant in the 
immediate neighborhood of the grazing orbit. 


4.2 Investigating the Trace for Singularity 

To obtain the expression for the trace of the Jacobian J in d?]), we need to obtain first the expressions for 
the maps Pi and P 2 - 

Let us consider a periodic solution to the equation given in (|9]) as 

p{t) = {u{t),v{t)f. 


Let {p(t) + Sp(t)) be a perturbed orbit, where Su satisfies the following variational equation 

Sil + 2(uJnSu + OJn^Su = 0. (13) 

The variational equation needs to be solved to obtain the perturbed flow 6 p{t) = (dit(r), du(r))^. Solving 
the variational equation amounts to solving the first-order differential equations 

d ( Su{t) 
dt I Svlt) 


0 1 


6u{t) 

Svlt) 


with (5u(0) = Suo,Sv(0) = Svq. 

The solution of the above problem can be expressed as 


where 


Nr = e 


— „-C‘^nT 


cos{u:qt) 


Su{t) 

SvIt) 




= Nr 


sin(wo'r) 


Suo 

Svq 


sin(woT)/a;o 




--U!q sin(a;o'r) cos(wor) — 




: sin(a;oT) 


(14) 


with Wo = U!n\/l — C^- 

Now we can proceed to obtain the expression of the trace of the Jacobian in ([7]). In the situation shown 
in Fig. 01 the flow takes time sq to reach the discontinuity boundary starting from the Poincare plane, and 
time {T — sq) to return to the Poincare plane starting from the discontinuity boundary, where T is the time 
period of the external forcing function g(t). Thus using the notation in (I14L 


Ni 

KJ' 

1 

II 

0 

II 

N 2 

= N(^t-so) = 


nil rii2 

ni3 ni4 

n2i n22 
n23 n2A 

where the expressions for nn, ni2, ni3, ni4, n2i, n22, ^23 and n24 can be derived from ( 03 . 
Using ( 1121 ), (Uni) and (HH]) in (I?]), we get 

10 
a 1 


(15) 

(16) 


J = 


n2i n22 

n23 ^24 


nil ni2 
ni3 ni4 






n2i J- an22 n22 
n23 + Q!n24 n24 


nil ni 2 
ni3 ni4 


^211^,11 + n22ni3 -|- an22nii 


« 23«12 + n24ni4 + an24ni2 

Tr{J) = e“‘’‘^"^{n 2 inii-I-n22ni3 J-n23ni2 + n24ni4 J-a 


{n22nii + n24ni2) } 


(17) 










The expression for the trace of the Jacobian Tr{J), in (I17L shows that the singularity term a has a 
coefficient e“'“‘^"^(n 22 nii + n 24 ni 2 )- Let us take a closer look at this coefficient. Using (ITTl) . (ITS|) and (fTBl) 


nil = cos(a;oso) + 


c 




: sin(a;oSo) 


ni2 = 


«22 = 


sin(a;oSo) 

Wo 

sin{a;o(r- sp)} 
Wo 


and n24 = cos{wo(r — so)} — 


x/T^ 


: sin{wo(L’ - So)}. 


Therefore it follows that 


I sin(a;oT’) / rx 

^ 22^11 + n2Ani2 = —^ 0 


V Wo 7^ 


'^Wforcing 

2 


(18) 


where Wforcing is the angular frequency of the periodic forcing function g(t), i.e., WforcingiL = 27r, and m € I. 
Thus the coeffcient of a in the expression of the trace (II3 of the Jacobian of the stroboscopic map must be 
a non-zero entity. Thus the singularity in a survives, and hence a square-root singularity must occur in the 
trace of the Jacobian. 


5 Conclusions 

In this paper we have probed the variation of the trace and determinant of the Jacobian matrix of map of a 
hard-impact oscillator as it goes from non-impacting state to an impacting state. We have shown that the 
square-root singularity should be expressed only in the trace of the Jacobian matrix while the determinant 
should remain invariant in the immediate neighborhood of a grazing orbit. 
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